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Abstract. Polarization measurements for the optical counterpart to GRB 021004 are presented and discussed. Our 
observations were performed with the TNG and the VLT-UT3 (Melipal) during the first and fourth night after the 
gamma-ray burst discovery. We find robust evidence of temporal evolution of the polarization, which is therefore, 
at least partially, intrinsic to the optical transient. We do not find convincing evidence of wavelength dependence 
for the intrinsic polarization of the transient, in agreement with current polarization models for optical afterglows. 
We discuss the role of dust, both in our galaxy and in the host, in modifying the transmitted polarization vector, 
showing how a sizable fraction of the observed polarized flux is due to Galactic selective extinction, while it is 
not possible to single out any clear contribution from dust in the host galaxy. We discuss how our data compare 
to those obtained by different groups showing that a two-component model is required to describe the complete 
dataset. This is not surprising given the complex lightcurve of GRB 021004. 
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1. Introduction rule out that the observed polarization was induced^ by 
^ dust along the line of sight in our own Galaxy (e.g. Covino 
- - It is now well established that gamma-ray burst (GRB) et al. 119991 r 2UU3c|l . More recently, it was also possible to 
optical afterglows (OA) can show some degree of linear exclude, at least for some cases, a major contribution to 
polarization. To date, in five cases a positive detection was the observed polarization level due to the interposition of 
obtained: GRB 990510 (Covino et al. 119991 Wijers et al. dust in the host galaxy. In fact, interstellar polarization is 
UnSIl, GRB 990712 (Rol et al. lMTUjl . GRB 020405 (Bersier necessarily associated with reddening, and by modeUing 
et al. |2D03a, Covino et al. I2003al Masetti et al. J003), the spectral shape of the OA, the total amount of dust in- 
GRB 020813 (Barth et al. 120031 Covino et al. I2002a|l . terposed on the line of sight can be constrained, allowing 
and GRB 030329 (Covino et al. I2U03I Efimov et al. 120031 us to put limits on the dust-induced polarization. 
Magalhaes et al. 2003). Usually, the polarized flux is not This method is however model dependent, and re- 
large, P <, 3%, but in most cases it has been possible to lies also on the knowledge of the dust properties in the 

GRB environment (see e.g. Lazzati et al. I2002a)l . Given 
the present uncertainties, the possibility that for at least 

'7, Z a ~^ 4- J- T\ T 4-- -1 some OA a sizable fraction of the observed polarization 
bend offprint requests to: D. Lazzati; e-mail: 

lazzati@ast cam ac uk induced by dust in the GRB environment or in the 

* Based on observations made with ESO telescopes at the 

Paranal Observatory under programme Id 70.D-0111, on data ^ Here and in the following the only dust induced polariza- 

from the ESO/ST-ECF Science Archive Facility and on obser- tion that we consider is the one due to the dichroism of the 

vations made with the TNG under programme TAC8_01(47). aligned grains and not the polarization induced by scattering. 
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host galaxy cannot be excluded yet. In principle, there 
are at least two safe ways to unambiguously detect in- 
trinsic polarization: the first is to perform multiple obser- 
vations, looking for temporal variation of the polarization 
degree and/or position angle. Such a polarization variabil- 
ity would also provide a direct link between the dynamics 
of the fireball evolution and the geometry of the emitting 
region (Ghisellini & Lazzati I1999L hereafter GL99; Sari 
119991 Granot et al. l2()02(l . The second is to study the wave- 
length dependence of the polarization, possibly extending 
it to the infrared, in an attempt to exclude that it follows 
the "Serkowski curve" typical of Milky Way (MW) inter- 
stellar dust polarization (Serkowski et al. I1975|l . bearing 
in mind, however, that the OA emission could be intrinsi- 
cally polarized in a wavelength dependent way. 

Recently, some steps in this direction were performed. 
First, for GRB 020813 secure polarization variability was 
detected, the degree of polarization decreasing on a day 
time scale from P ~ 2% (Earth et al. .2003,1 to P ~ 
0.8% (Covino et al. I2002a|l . at a fixed polarization angle. 
Moreover, again for GRB 020813 (Barth et al.EMI, and 
very recently for GRB 030329 (Covino et al. lM75|l . spec- 
tropolarimetry could be performed; in both cases, small 
but significant wavelength dependence was found. 

GRB 021004 was localized on 2002 October 4 at 
12:06:14 UT by the HETE~II satellite (Shhasaki et al. 
Onn^ . In the FREGATE 8-40 keV and in the WXM 2- 
25 keV bands the burst had a duration of about 100 sec- 
onds. It thus belonged to the class of long-duration GRBs. 
The optical counterpart was identified less than 10 min 
after the burst (Fox I20()2|l as an P ~ 15.3 fading ob- 
ject at the coordinates a2ooo = 00''26™54!69, (52ooo — 
-hl8°55'41'.'3. The early detection of the OA and its bright- 
ness allowed a dense sampling of the light curve (see e.g. 
Lazzati et al. I2002bl and references therein) , the identifi- 
cation of several absorption systems, and of a prominent 
emission feature in the optical spectra identified as a Lya 
emission line from the host galaxy at a redshift z = 2.328 
(Mirabal et aLEOOS ?; Matheson et al.EOnS M0ller et al. 
Schaefer et al. 11003) • 

In the following we will discuss our three polarimetric 
observations of GRB 021004: one performed in the near in- 
frared (NIR) with the Italian Telescopic Nazionale Galileo 
(TNG) at the Canary Islands, and two in the visible band 
at the VLT. Furthermore, we analyzed a publicly available 
spectropolarimetric dataset retrieved from the ESO VLT 
Science Archive. We will then compare these measure- 
ments, including the polarimetric observations performed 
by Rol et al. ((2003(1 . with theoretical models. 




0.01 0.1 1 10 100 

Time since GRB (days) 

Fig. 1. P-band hghtcurve of GRB 021004 with the posi- 
tion of the polarimetric observations marked. Data are 
from Bersier et al H2003a;i : Fox et al. ijMT^ : Holland et 
al. iPMjl : Pandey et al. 1^70^ : Uemura et al. l|M75j> . 



imaging polarimetry mode, and with the grism 300 V in 
the spectropolarimetry mode. 



2.1. TNG observation 

The TNG observation (hereafter run 1) started on October 
4.915 (9.86 hours after the GRB trigger) and lasted for 
^ 1.8 hours. The optical transient (OT) J magnitude was 
derived by the acquisition frames as J = 17.00 ± 0.05 
(hereafter l-cr errors are reported) . The observations were 
performed under mediocre seeing conditions (1.5") in the 
large field mode with a scale of 0.25"/pixel. 



2.2. VLT observations 

We analyzed three sets of polarimetric observations of 
GRB 021004 carried out by the ESO VLT-UT3 (Melipal). 



2. Observations and analysis 

Observations of GRB 021004 made use of two tele- 
scopes. First, the TNG, equipped with the Near Infrared 
Camera Spectrometer (NICS) and a J filter, in the 
imaging polarimetry mode; second, the ESO VLT-UT3 
(Melipal), equipped with the Focal Reducer/low disper- 
sion Spectrometer (FORS 1) with a Bessel V filter in the 



First VLT observation 

Our first VLT observation (here called run 2) started on 
October 5.080 (13.82 hours after the GRB trigger) and 
lasted for ~ 1.6 hours (see Fig.^. The OT was clearly 
detected in the acquisition image with a magnitude V — 
19.34 ±0.02 with respect to the USNO-A2.0 star reported 
by Fox <|^nn^ . as calibrated by Henden ((2002a. ■2002b|l . 
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Second VLT observation 

From the ESO archive we retrieved a pubhc spectropolari- 
metric observation (run 3). It started on October 5.247 
(17.83 hours after the GRB trigger) and lasted for ~ 2.0 
hours. The spectrum covers the range from 350 nm to 
860 nm and the observations were performed with the 
300 V grism. 

Third VLT observation 

Eventually, a last VLT observation (run 4) was performed 
starting on October 8.225 (89.3 hours after the GRB trig- 
ger), and lasting for 2.8 hours. The OT magnitude was 
V = 20.89 ± 0.03 with respect to the same USNO-A2.0 
star. 
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Fig. 2. Polarization degree, position angle and normalized 
spectrum of the public ESO-VLT spectropolarimetric ob- 
servation (run 3). The polarization degree is roughly con- 
stant throughout the optical band at P ~ 1.9% and the 
position angle is also constant at ~ 118°. Observations 
were performed with the VLT-UT3 (Melipal) equipped 
with FORS 1 in the spectropolarimetry mode with grism 
300 V. Filled symbols correspond to the "clean" part of the 
spectrum, while empty symbols are relative to absorbed 
portions of the spectrum, or the ones contaminated by 
strong sky emission lines. 



All VLT observations were performed under 
good/ excellent seeing conditions (0.5" — 0.9") in standard 
resolution mode with a scale of 0.2"/pixel. 

Polarimetric standard stars were also observed. One 
polarized, BD-125133, in order to fix the offset be- 
tween the polarization and the instrumental angles, and 
three non-polarized, WD 0310-688, WD 1615-154, and 
BD-l-284211, to estimate the degree of spurious polariza- 
tion possibly introduced by the instruments. In addition, 
we have also analyzed the ESO archive spectropolarimen- 
tric data for the NGC 2024 NIRl polarization standard 
star, obtaining in all cases a good Serkowski curve com- 
pletely consistent with the available data^. 

2.3. Data reduction and analysis 

The data reduction was carried out with the Eclipse 
package (version 4.3.1; Devillard I1997|l . After bias sub- 
traction, non-uniformities were corrected using flat-fields 
obtained without the WoUaston prism. For the IR data 
sky flat-fields were applied. The flux of each point source 
in the field of view was derived by means of aperture pho- 
tometry by the Graphical Astronomy and Image Analysis 
(GAIA) tools^ (version 2.6-6). 

The general procedure followed for the analysis of 
imaging polarimetry observations is indeed extensively 
discussed in Covino et al. I2002bl I2()()3ajl and di 

Serego Alighieri fl997!), while details about the NICS po- 
larimetric capabilities are discussed by Oliva (1997). 

The VLT Spectra were extracted with the apall tool 
included in the IRAF package (version 2.11), which ex- 
tracts the spectrum in a fixed-width window, allowing for 
a polynomial evolution of the centroid with wavelength. 
Suitable IDL routines were also developed to perform an 
independent extraction and check for any possible bias. In 
this case the extraction of the ID spectrum from the 2D 
frame was performed by fitting a bell-shaped function to 
all the vertical stripes of the 2D spectrum, allowing for a 
first order variation of the background and for a 4*^ de- 
gree polynomial variation of the flux centroid and width of 
the fitting function. The result of this "fitted" extraction 
were not entirely consistent with those of a more standard 
sliding window, as discussed below. 

The results of the spectropolarimetric observations are 
shown in Fig.|21 To all our data we have applied a cor- 
rection for the wavelength dependence of the half-wave 
plate axis as measured by the FORS team*. The plot is 
based on the fitting extraction (hereinafter F) which, as 
anticipated, is not fully consistent with the more stan- 
dard window extraction (hereinafter W). The F extrac- 
tion shown in the figure yields a polarization which is 
consistent with being wavelength independent. A simul- 
taneous constant fit to the polarization degree and an- 
gle returns P = (1.88 ± 0.05)% and ■& = (118 ± 1)° 



^ http : //www. eso . org/ instrument s/f or si /pola. html 

^ http : / / star-www . dur . ac .uk/ ~pdraper/gaia/gaia.html 

^ http : //www. eso . org/ instrument s/f orsl/ 
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with xVd-O.f. = 45.8/44 (null probability ~ 60%). The 
bin size has been made large enough to make the po- 
larization bias (Wardle & Kronberg |1974| unimportant. 
Completely consistent results are obtained with the W 
extraction if a small window (6 pixels width) is adopted, 
while a wider window (20 pixels) yields different results. 
The polarization obtained using the wide W extraction is 
wavelength dependent, the polarization decreasing from 
~ 2% at 600 nm to ~ 1% at 800 nm. The resulting polar- 
ization also presents a significant dip at P ~ 1% around 
550 nm. In principle, both methods are prone to inadequa- 
cies in extracting the data, and it is not possible to decide 
a priori which is the best one. The fit and small W extrac- 
tion are in fact sensitive to non centrally symmetric 2D 
spectra. The large window, instead, can yield an unbiased 
extraction of non-symmetric 2D spectra, but is subject 
to a varying background and to contribution from nearby 
sources. In an "Ockham razor" approach, we decided to 
show the result which requires less parameters to be ex- 
plained, namely the constant one, warning however the 
reader that according to the large window extraction, the 
polarization may be wavelength dependent, being smaller 
in the red part of the spectrum. 

From an independent analysis of the same data, Wang 
et al. H2()()8|l report a marginal evidence of increasing po- 
larization at A ^ 400 nm, across the rest-frame Lya ab- 
sorption features. They interpret this as a consequence of 
absorption of the fireball emission by nearby high veloc- 
ity clumps of HI atoms, with a covering factor smaller 
than, but comparable to, unity. As can be seen in the up- 
per panel of Fig. |21 our binned data do not support this 
claim, and the same result is obtained with a binning of 
~ 2 nm, comparable to the one adopted by Wang et al. 
PMjl . However we find that using smaller wavelength 
bins (e.g. 0.52 nm) the polarization, not corrected for the 
bias (Wardle & Kronberg I1974|l . steadily increases from 
400 nm down to the shortest available wavelength as ex- 
pected for a low sensitivity region at the blue edge of the 
spectral range. Therefore this UV rise of the polarization 
looks consistent with a bias effect, although the possibil- 
ity that some of it might be real cannot be completely 
excluded. 

We stress that the difference between the various spec- 
trum extraction methods discussed here are at the ~ 1% 
level. While being small, if not negligible, for a stan- 
dard spectral analysis, these differences become significant 
when weak polarization levels are investigated. 

2.4. Other observations 



Rol et al. H2()();-{|l obtained three independent measure- 
ments of the optical polarization of the afterglow of 
GRB 021004, three with the Nordic Optical Telescope, 
in the i?-band, and one with the VLT, in the l^-band. 
They obtain P = (1.17 ± 0.46)% with position angle 
I? = (184.2 ± 11.4)° at At = 0.366 d after the GRB, 
P = (1.73 ±0.51)% with position angle -d = (166.4 ±8.1)° 
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Fig. 3. V^-band time-resolved polarimetry of GRB 021004. 
Filled circles show the polarization degree (upper panel) 
and position angle (lower panel) for our VLT runs (run 
2 and run 4) vs. the time elapsed since the GRB trigger. 
The diamonds refer to the Rol et al. H2003|l measurements, 
while the asterisk is obtained by filtering with a Bessel 
V transmission profile the spectropolarimetric observation 
(run 3). The shaded bands show the MW interstellar po- 
larization for field stars. Their vertical width corresponds 
to the l-cr uncertainty. 



at At = 0.376 d, P < 1% at At = 0.396 d and P = 
(1.29 ± 0.13)% with position angle 'd = (121.8 ± 2.8)° at 
At = 0.666 d. 



3. Results and modelling 

3.1. Imaging polarimetry 

The IR polarimetric observation (run 1) provided Q and U 
Stokes parameters compatible with those derived for run 2 
(see below and Tab.^, even though the larger errors only 
allowed us to derive a P < 5% upper limit (95% confi- 
dence level). This limit is not particularly stringent com- 
pared to those derived in the optical for other afterglows 
(e.g. Hjorth et al. .l999. Covino et al. I2002bl Bjornsson et 
al.Hnnil Covino et al. l2003d|l . However, it is the first IR 
polarimetric observation that provides a useful constraint, 
since the previous upper limits were rather loose (~ 30%; 
Klose et al.|MD- 

The results of VLT l^-band imaging polarimetry, as 
well as a synthetic 1^-band measurement derived from the 
spectropolarimetric measurement, are plotted in Fig. |31 
as a function of the time elapsed since the GRB trigger. 
We also report the value of polarization induced by the 
MW selective extinction, as derived by averaging the Q 
and U parameters for several bright stars in the field. We 
obtain Pmw = (0.81 ± 0.03)% and i9mw = (107 ± 1)°. 
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Table 1. The normalized polarization Stokes parameters not corrected for the interstellar polarization. Observations 
were performed with the TNG (run 1) and with the VLT-UT3 (runs 2 and 4). VLT-UT3 spectropolarimetry (run 3) 
was performed with grism 300 V. The reported results for the spectropolarimetric observation of run 3 were obtained 
by integrating the spectrum over the V^-band. Uncertainties are at 1-a and the upper limit is a 95% confidence level. 



Run 


Tel. 


Filter 


UT (10/02) 


Magnitude 


Q 


u 


P (%) 


^ n 


1 


TNG 


J 


4.953 


17.00 ±0.05 


-0.0160 ±0.0130 


-0.0222 ±0.0130 


< 5 




2 


VLT 


V 


5.172 


19.34 ±0.02 


-0.0083 ± 0.0009 


-0.0094 ± 0.0010 


1.26 ±0.10 


114 ±2 


3 


VLT 


V 


5.247 




-0.0075 ± 0.001 


-0.016 ±0.001 


1.74 ±0.20 


122 ±2 


4 


VLT 


V 


8.225 


20.89 ± 0.03 


-0.0067 ± 0.002 


±0.0002 ± 0.002 


0.67 ±0.23 


89± 10 



Such a polarization is slightly larger than the maximum 
interstellar (ISM) polarization according to the empirical 
law PisM ^ 0.09 Eb-v (Serkowski et al. [T775|l . given the 
reddening Eb-v ^ 0.06 estimated for the GRB 021004 
field (Schlegel et al. I1998|l . We note however that devi- 
ations from the above law are observed especially along 
low-reddening lines of sight (see Fig. 9 of Serkowski et al. 
I1975|l . An important conclusion that can be drawn from 
Fig. |3| is that neither the polarized fraction nor the posi- 
tion angle of the OT were constant during the evolution of 
the afterglow (Rol et al. I2003|l . We conclude that the rel- 
ative contribution of the polarizing components, namely 
the MW ISM, the GRB host ISM, and the OT itself, did 
vary on short time scales. Since the ISM polarization can- 
not vary on short (day) time scales, we conclude that the 
OT was intrinsically polarized. The rotation of the polar- 
ization angle between the first measurements of Rol et al. 
H2(JU3|I and our late-time datum is consistent with 90°, a 
quantity predicted in the more commonly accepted mod- 
els for the polarizations of GRB jets (GL99; Sari 11999(1 . 
However, in contrast with the prediction, the transition 
between the two angles is not sharp, but rather smooth. 
This suggest the presence of an additional polarizing com- 
ponent, as we will discuss below. 

3.2. Spectropolarimetry 

We have shown in Fig. |3 the result of our reduction of 
the spectropolarimetric observation publicly available at 
ESO. Spectropolarimetric observations of OTs are impor- 
tant since they allow us, in principle, to single out which 
component is contributing more to the observed polariza- 
tion: the OT itself, the host ISM, or the MW ISM. This is 
possible since the three components have different wave- 
length dependencies. 

The OT polarization is supposed to be wavelength in- 
dependent, at least in the limited spectral range investi- 
gated here. ISM polarization is instead wavelength depen- 
dent. It follows a "Serkowski law" : 



PW = Pmax exp 



.K In^ 



(1) 



where Pmax is the maximum induced polarization, ob- 
tained at Amax- Experimentally, 0.34 /im ^ Amax ^ 
0.9 fim, and this parameter is considered to be a mea- 
sure of the size of the polarizing grains: the larger Amax, 



Table 2. Results of the simple modelling of the spectropo- 
larimetric result. ^ The constant polarization value is given 
for the OT model. ^ The value Amax = 0.34 ^m has been 
set as a lower limit to the parameter (see text). 



Model 


Pmax(%) (°) An,ax (/im) x7d.0.f. 


MW ISM 
Host ISM 


1.88 ±0.05 118 ±1 — 45.8/44 
1.95 ±0.1 118 ±1 0.52 ±0.05 50.5/43 
2.1 ±0.15 118 ± 1 0.34* 57.2/43 



the larger the grains. For the numerical value of the coef- 
ficient K, we follow the more recent study of Martin et al. 
l(T^ : 



K = 



1.66 Amax 

-0.59 ±2. 56 A„ 



A > A„ 
A < An 



(2) 



To model the polarization results shown in Fig. |21 we 
have assumed first that the entire polarization is due to 
only one of the three possible components (i.e. OT, MW, 
host), taking into account the wavelength dependence of 
the ISM polarization as detailed above, and calculating the 
effect of the host-ISM in the host rest frame^. As reported 
in Tab. 12 the quality of the data allows us to exclude a 
dominant role of the host ISM in producing the observed 
polarization, but it is not possible to identify a dominant 
component between MW-ISM and intrinsic. The first two 
fits are both acceptable (x^/d.o.f. < 1.2), with a slight 
preference for the OT model, while the host-ISM fit yields 
a significantly worst x^- Performing the fit only on the 
high quality dataset (filled symbols in Fig. |2J) yields sta- 
tistically indistinguishable results. However, since the field 
stars do show a moderate degree of polarization, we know 
that the observed polarization cannot be entirely due to 
the OT itself. Furthermore, at all frequencies the polar- 
ization position angle is very similar to that of the field 
stars in the imaging polarimetry (see above). 



^ The application of the Serkowski law to high-redshift sys- 
tems is not directly supported by observations, given the in- 
trinsic difficulty in performing spectropolarimetric studies of 
distant objects. It is however at least partly justified by the 
observation of interstellar polarization following a Serkowski 
curve in SN 1986G in Centaurus A, at z = 0.00183 (Hough et 

ai. nwt . 
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Fig. 4. Effect of a polarizing ISM on the theoretical polarization of afterglows. The theoretical curves (top: polarization 
level; bottom: position angle) from GL99 for 'do/'&j — 0.5 are shown, in both panels, as solid lines (pisM = 0). Lines with 
different styles show the effect of the selective extinction for ISM with small, intermediate and comparable polarization 
with respect to the OT maximum polarization. The position angles of the ISM are reported relative to the initial angle 
of the OT. The left panel shows the effect of an ISM that induces polarization with a position angle i?ism = —19°, 
while the right panel shows an ISM with i?ism — 60°. 



3.3. Combined modelling 

Since the polarization levels are small and different effects 
seem to contribute at comparable degree, it is mandatory 
to perform a combined modelling of the imaging and spec- 
tral polarimetry, combining the effects of the OT polariza- 
tion and of the MW ISM selective extinction. The need of 
at least two polarizing components is due to the presence 
of time variability of the measured polarization (OT 
component) and b) of significant polarization in the field 
stars (MW ISM component). 



3.3.1. Transmission of variably polarized light in a 
polarizing ISM 

In order to combine the effects of ISM selective extinc- 
tion with the evolving intrinsic polarization of the OT, 
we adopt a Mueller calculus approach (see e.g. di Serego 
Alighieri 119971 and references therein). In this formalism, 
the transmitted Stokes vector S' = (/', Q' ,U',V') is com- 
puted from the incident one S = (/, Q, C/, V) through a 
matrix, called "Mueller matrix" , which incorporates all 
the properties of the transmitting medium: S' = M ■ S. 
Since developing a complete treatment of the polarizing 
properties of the ISM is far beyond the scope of this paper, 
here and in the following we adopt a simplified version of 
the Mueller matrix which is correct for a non-birefringent 
dichroic medium. With this simplification we assume that 
the ISM is not able to induce circular polarization to any 
incident light, either polarized or not. In fact, the ISM 



does induce a small degree of circular polarization even in 
unpolarized sources (Martin & Angel [T976|l . which indi- 
cates that a non-coaxial birefringent and dichroic medium 
should be considered. However, the optical properties of 
such a medium are described by five parameters (extinc- 
tion, orientation of the optical axes of birefringence and of 
dichroism and the two respective refraction indexes) while 
we can constrain only four parameters observationally (ex- 
tinction, induced linear polarization, position angle and 
induced circular polarization). It is therefore not feasible 
to derive a complete Mueller matrix without going into a 
detailed modelling of the structure and geometry of dust 
grain and their alignment. 

Consider a dichroic ISM that induces a polarization 
piSM = {q^ + it^)^/^ on unpolarized stars. Its Mueller ma- 
trix has then the form: 



M = 



/I 
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PiSM 


u 


qu{\-A) 
PlSM 


\o 






u 
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where A = -^/l — Pigj^, and e^'^ is the opacity of the 
medium to non-polarized radiation. In order to conserve 
energy (i.e. not to have an increased transmitted inten- 
sity), this parameter must satisfy: 

1 



< 



1 



(4) 



the equality holding for a perfect polarizing medium, i.e. 
one that does not absorb any radiation completely po- 
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larized with its same angle. It can also be easily shown 
that if pisM ^ 1 and Q,U,V ^ 1 the matrix com- 
putation is equivalent to a simple sum of the incident 
and ISM normalized Stokes parameters: Q' « q + Q, 
U' ~ u + U. Formally, the Mueller matrix in Eq. O is 
not valid in the case of a birefringent medium. It can 
however be shown that it is a good approximation (de- 
viations < 10%) if the incoming light is moderately po- 
larized {P < 10%), the linear and circular interstellar po- 
larizations are small (pisM < 10% and wism/pism < lO^'^; 
Martin & Angel I1976|l and the circular polarization in- 
duced on polarized sources is small {V/P < 10~^ for ac- 
tive galactic nuclei; Landstreet & Angel [T^72(l . 

In order to exemplify the effect of a polarizing ISM on 
the intrinsic OT polarization, we adopt the polarization 
model by GL99. This model predicts a polarization curve 
characterised by two distinct peaks, whose absolute inten- 
sity depend on the angle between the line of sight and the 
axis of the fireball (assumed to be coUimated in a jet). 
The polarization position angle is rotated by 90 degrees 
at the time of null polarization between the two peaks. 
This model is illustrated by the solid lines in Fig. ^ (i.e. 
no ISM polarization). With different line-styles, this figure 
also shows how the predicted polarization degree and po- 
larization angle change in time once modified by some in- 
tervening polarizing ISM. All angles are reported relative 
to the initial OT polarization angle. In the left panel the 
case relative to this burst is shown = —18°; see below), 
while in the right panel a larger misalignment between the 
OT and ISM angles is shown. The figure shows that the 
presence of the ISM modifies quite substantially the ob- 
served polarization, especially at times when Pqt ~ ^'ism- 
Note also that the presence of the polarizing ISM makes 
the position angle to vary smoothly, instead of sharply. 





i^— , — , — , — .— — 1 , , , — 


/ J 

/ : 
/ : 

/ : 
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Fig. 5. Time resolved modelling of the polarization and 
position angle of GRB 021004. The first three and the fifth 
datum are taken from Rol et al. H2003fl . Our NIR upper 
limit is not included. The solid line show the best fit ob- 
tained from the whole dataset, with a free break time and 
free properties of the ISM polarization. Even with this 
extra freedom a good fit cannot be obtained due to the 
rapidity of the evolution of the polarization (see text for 
more details). The dashed line shows instead the best fit 
model obtained by modelling only the last four data. In 
this case the model yields acceptable x^/d-O.f. — 6.3/4 
and a best-fit break time comparable to the one derived 
from the lightcurve fitting (see text for more details). The 
vertical dotted line shows the break time tb = 4.74 d de- 
rived from lightcurve modelling (Holland et al. 12003(1 . 



3.3.2. Modelling of the polarization curve 

To model the time-resolved polarization measurements in 
the framework of available models (GL99) with a polariz- 
ing ISM we perform a fit to the polarization curve propa- 
gating the intrinsic OT polarization through the MW ISM 
with the use of the Mueller matrix derived in Eq. |3| The 
model has 6 degrees of freedom; the initial position an- 
gle of the OT intrinsic polarization, the off-axis angle of 
the line of sight, the degree of alignment of the magnetic 
field, the jet break time and the ISM q and u parameters. 
The last three parameters can be constrained with the ob- 
servations, but we let them free to vary here to check a 
posteriori the agreement of the derived values with the ob- 
servations. In addition, we adopt models with and without 
sideways expansion of the jet (GL99; Sari ll99(?l Rossi et al. 
2003, in prep). A particularly important issue is to check 
whether the prediction of the models are consistent with 
the position angle rotation detected by Rol et al. H2003|l . 
The fit is performed on the q and u parameters rather 
than on p and due to their better statistical properties. 



The result of the fit is that it is not possible to model 
the ^ 90° rotation of the position angle in the framework 
of the proposed models, not even with the addition of a 
polarizing ISM with free properties. In fact, the smallest 

that can be obtained is — 43 for 8 degrees of free- 
dom (for a non sideways expanding jet). In addition, the 
fit formally yields a best break time tj = 0.25 d, in dis- 
agreement with the value tj — 4.74lQ g'* inferred from the 
lightcurve modelling (Holland et al. '20031. Also, the best 
fit excludes the presence of a sizably polarizing ISM, in 
disagreement with the observation of polarization of the 
stars in the field of GRB 021004. This could be explained 
with a polarizing ISM in the host galaxy with properties 
opposite to those of our ISM. Such a possibility would 
however require an unacceptable fine tuning and is ruled 
out by the modelling of the spectropolarimetric observa- 
tion (see § 3.2). The formal best fit model, re-converted in 
polarization and position angle, is shown in Fig. [S] with a 
solid line. 

There is one important issue, however, that makes this 
result non conclusive. The models for polarization that we 
have used in the fit above are derived under the assump- 
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tion that the firebaU and the ISM are homogeneous. In 
this case, the total lightcurve should be characterised by 
a smooth broken power-law decay. This is not the case for 
GRB 021004. Its lightcurve, as shown in Fig. ^ (see also 
Lazzati et al. l2002b|l . shows prominent bumps overlaid on 
this power-law. There are several possible explanations for 
these bumps: inhomogeneities in the ISM (Lazzati et al. 
I2002b|l . inhomogeneities in the fireball (Nakar et al. l2003|l 
or delayed injection of energy from the central engine (Fox 
et al. l2()();-i(l . While in the latter hypothesis the polariza- 
tion curve should not be affected, inhomogeneities, either 
in the fireball or in the ISM, can produce a polarization 
signal by breaking the symmetry of the fireball emission 
at early times (Granot & Konigl 2003). Such polarization 
would have a randomly oriented position angle, depend- 
ing on the azimuthal location of the inhomogeneity with 
respect to the line of sight. A similar effect is predicted in 
the case of micro-lensing events (Loeb & Perna [TI)98|) . If 
we look more closely at Fig. we note that the first three 
measurement of Rol et al. (|2()()3|l lie on top of a major 
rebrightening, while the remaining four observations were 
performed in relatively unaffected time intervals. Only the 
last four measurements should therefore closely follow the 
theoretical model. 

To check this idea, we have performed the same mod- 
elling described above only on the final four points, fixing 
the ISM polarizing properties to that of the stars in the 
field of GRB 021004 in order to limit the number of free 
parameters. Indeed, the last four points can be successfully 
described by the model, with a break time tj = 3± 1 days, 
in good agreement with the tj = 4.74^Q g'* d obtained by 
Holland et al. P00 3) from the break in the lightcurve. The 
best fit model, which has = 6.3 for 4 degrees of free- 
dom, is shown with a dashed line in Fig.jS] The fit requires 
a large degree of alignment of the magnetic field of 75% 
and a moderate off-axis line of sight i^o/^^j = 0.45. The 
best fit model is obtained for a non sideways expanding 
jet. A sideways expanding jet cannot however be rejected, 
with x^/d.o.f. = 7.2/4. This ambiguity is due to the fact 
that, according to the best fit model, the measured points 
lie in the first peak of the polarization curve, where the 
sideways expansion has only a marginal effect. The best fit 
model requires a misalignment of (—19 ± 1)° between the 
OT initial position angle and the ISM polarization. The 
effect of such a misalignment is shown in the left panel of 
Fig.H 

3.4. Complications 

Additional complications can be envisaged in locally red- 
dened afterglows. If in fact a strong polarization is induced 
by dust in the close vicinity of the burst explosion site, 
two effects can take place. First, dust destruction by the 
burst prompt and afterglow emission (Waxman & Draine 
120001 Perna & Lazzati [2002(1 can imprint a strong tempo- 
ral evolution in the host ISM polarization. This is however 
unlikely to be relevant for our observations, since the typ- 



ical time scale for this process is of the order of minutes, 
rather than days. Second, a patchy absorber may obscure 
in a different way different portions of the fireball (e.g. 
Wang et al. l2003|l . altering one of the key assumptions of 
the models. None of these effects should alter dramatically 
the lightcurve, but a random noise could be overlaid on the 
smooth theoretical evolution of the polarization predicted 
by the models. 

In addition, a preferential direction for the polariza- 
tion, different from the one defined by the plane contain- 
ing the jet axis and the line of sight, can be defined by 
the presence of an interstellar magnetic field of sufficient 
magnitude. Such a possibility has been recently studied by 
Granot & Konigl H2003|l . Qualitatively, the effect of such 
a pre-existing field is not different from that of a polar- 
izing ISM. Due to the presence of a second asymmetry, 
the evolution of the polarization angle is smooth rather 
than sharp, and the external component dominates the 
observed polarization properties when the intrinsic OT 
one is small. Such additional degrees of freedom were not 
necessary in our last fit (nor they can explain the early 
~ 90° angle rotation). We stress however that, since all 
our data were observed before the jet break, we cannot 
exclude the presence of an ISM magnetic field with a well 
defined orientation. 

4. Summary and conclusions 

We have presented multi-time and multi-filter observa- 
tions of the polarization of the afterglow of GRB 021004 
performed with ESO-VLT and TNG as well as our anal- 
ysis of the publicly available ESO-VLT spectropolarimet- 
ric observation. The interpretation of the observations is 
complex since none of the polarizing mechanisms that con- 
tribute to the observed polarization seems to clearly dom- 
inate over the others. We therefore adopt this afterglow as 
a case study to investigate and describe these three main 
effects: intrinsic (and time varying) OT polarization, host 
ISM polarization, and MW ISM polarization. By mod- 
elling the spectropolarimetric and time-resolved imaging 
polarimetry we were able to get rid of the contribution of 
the host ISM, while OT and MW ISM polarizations seem 
to play an intertwined role, one dominating over the other 
at different times. 

To perform a detailed time dependent modelling of 
the polarization and position angle evolutions, we imple- 
mented our data-set with the four observations of Rol et 
al. ("2003). The complete dataset is particularly interesting 
since a sizable rotation of the position angle is present. Our 
attempt to fit the polarization data within the framework 
of GL99 models was however a failure, given the short 
timescale of the evolution and the lack of an appreciable 
break in the lightcurve at the time of the position angle 
evolution. Since the angle rotation is associated to one 
of the rebrightening events in the lightcurve, this burst 
is not suited for a comparison with models, which are 
computed for homogeneous fireballs producing featureless 
lightcurves. A possible explanation within exsisting mod- 
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els is that the early time polarization is dominated by a 
local enhancement of the fireball emission, which breaks 
the symmetry of the fireball producing a polarized signal. 
Indeed, the late time polarization data can be successfully 
and consistently described by the models. We therefore 
suggest that the lightcurve bumps are due to local events 
on the fireball surface, such as inhomogeneities in the ISM 
(Lazzati et al. |2002b ) or within the fireball itself (Nakar 
et al. 12003 ■ Refreshed shocks (Fox et al. '2003) are in- 
stead unable to account for the combined lightcurve and 
polarization evolution and can be rejected. 

The quality of the data, and in particular the lack of 
late time measurements) does not allow us to pin down 
the geometry and dynamic of the outflow, but we have 
shown that, with good quality spectropolarimetry and 
multi-time/multi- filter imaging polarimetry, it is in prin- 
ciple possible to disentangle the three effects and get a 
hold on the intrinsic polarization and on the structure 
and dynamics of GRB outflows. The added value of such 
a measurement would be the study of the polarizing prop- 
erties of dust in high redshift galaxies, a poorly studied 
property of such an important component of high redshift 
objects. 
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